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We consider a rotating neutron star with the presence of hyperons in its core, using an equation of 
state in an effective chiral model within the relativistic mean field approximation. We calculate the 
hyperonic bulk viscosity coefficient due to nonleptonic weak interactions. By estimating the damping 
timescales of the dissipative processes, we investigate its role in the suppression of gravitationally 
driven instabilities in the r-mode. We observe that r-mode instability remains very much significant 
for hyperon core temperature of around 10 s K, resulting in a comparatively larger instability window. 
We find that such instability can reduce the angular velocity of the rapidly rotating star considerably 
upto ~ 0.04 J1a', with Qk as the Keplerian angular velocity. 

PACS numbers: 97.60.Jd, 04.30. Db, 04.40.Dg, 26.60. +c 



I. INTRODUCTION 

The study of neutron stars are natural testing ground 
for studying extremely dense matter. The densities in 
the interior of such stars can reach up to several times 
the nuclear matter saturation density. At such high den- 
sities, with higher fcrmi momenta being available, high 
mass hadrons can be accommodated leading to a hyper- 
onic core in the neutron star interior p] |2] . There could 
also be the possibility that at such densities, when the 
nucleons are crushed, there could be quark matter [3J 
which can result in a color superconducting core in the 
interior of the neutron star [4j. In fact, there are dif- 
ferent possibilities of the ground state of dense matter, 
which could be stable strange quark matter [5] , and vari- 
ous possibilities of color superconducting matter [5] . The 
reason is that the true ground state of the dense matter 
system for densities relevant for the densities in the inte- 
rior of neutron stars is still an open problem because of 
the inherent nonperturbative nature of strong interaction 
physics. Moreover, external conditions like electrical and 
color charge neutrality conditions for the bulk matter in 
the interior of the star can also lead to various different 
possible phases of quark matter [3 [8] . This has given rise 
to various possibilities of compact stellar objects like neu- 
tron stars, strange stars, hyperonic stars or hybrid stars 
with a quark matter core and a crust of hadronic matter 
[2] • Obviously, it becomes very challenging to distinguish 
various compact stellar objects observationally. 

One of the various signatures suggested to distinguish 
different compact stars has been the r-mode instability 
[5HTTj. Various pulsating modes exist in neutron stars 
classified by the nature of the restoring forces. The r- 
modes correspond to the pulsating modes of the rotating 
stars where the restoring force is the Coriolis force and 
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these modes are axial modes [12]. As the fluid inside 
the star is self gravitating, these oscillations can couple 
to metric perturbations and lead to emission of grav- 
itational waves |13j . Gravitational modes drive the r- 
modes unstable due to Chandrasekhar-Fricdman- Schutz 
(CFS) mechanism [14j . This phenomenon will occur if 
the damping is sufficiently small and therefore provides 
a probe to study the viscosity of the matter in the interior 
of the star. While shear viscosity prevents differential ro- 
tation in a star, bulk viscosity dampens the density fluc- 
tuations in the star. Shear viscosity seems to be impor- 
tant at lower temperatures while bulk viscosity becomes 
the dominating dissipation mechanism at higher temper- 
atures (~ 10 9 A"). Further, since typical r-mode frequen- 
cies are of the order of rotational frequencies of the stars 
(Is -1 < to < 10 3 s _1 ), which are of the order of weak in- 
teractions, therefore the viscosity that is of relevance for 
r-mode instability will be dominated by the weak pro- 
cesses. Bulk viscosity is produced when the mode os- 
cillations induce perturbations in pressure and density 
and drives the system away from ^-equilibrium. As a 
result energy is dissipated from the system as the weak 
interaction tries to restore the equilibrium. While for 
hadronic neutron star matter, modified Urea processes 
(n+n — > n+p+e~ +u e ) involving leptons are important, 
it has been noted that the damping of the instability is 
dominated by large bulk viscosity arising due to nonlep- 
tonic processes involving hyperons [T5] . The correspond- 
ing viscosities are not only stronger, the temperature de- 
pendence is also different (varying as T~ 2 instead of T 6 
) which makes the hyperon bulk viscosity important at 
lower temperatures. Bulk viscosities of dense baryonic 
matter have been calculated under various assumptions 
as well as conditions over last few years [T51421] . 

The study of the r-mode oscillations provides an av- 
enue to study the density profile of the star along with 
the mass and the radius which in turn depends upon 
the equation of state for the matter in the interior of 
the star. Phenomenologically, for hadronic matter par- 
allel to the a — to model, popularly known as Walecka 
model [331 HI] i the chiral effective models have been de- 
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veloped and and are applied to nuclear matter. A chiral 
(7 — u> model along with dilatons in the context of phase 
transition has been developed in Ref. [25 and further gen- 
eralised to describe strange hadronic matter [2"(>H28|. An- 
other approach that was also considered is the parity dou- 
blet model [29] ■ We had recently considered attributes of 
rotating hyperonic star within a chiral a model |30H32j . 
This hadronic model also embodies a dynamical genera- 
tion of the vector meson mass along with nonlinear terms 
in the scalar field interactions to reproduce nuclear satu- 
ration properties at reasonable incompressibility [33J [M] . 
This model was then generalised to include the lowest 
lying octet of baryons [35 -37]. The resulting equation 
of state was used to calculate various gross properties 
of nonrotating as well as rotating neutron stars. It was 
further observed that the parameters of the model are 
sternly constrained and turn to be consistent with the 
constraints on nuclear equation of state from heavy ion 
collision data [3S1 EZ] • in the present work, we discuss 
the calculation of bulk viscosity within this model and its 
consequences regarding the r-mode instability problem in 
rotating neutron star with a hyperonic core. 

We organise the paper as follows. In the first sub- 
section of section II, we describe very briefly the chiral 
model that we shall be using to derive the equation of 
state needed to study the structure of the neutron star. 
In the subsequent subsections in the section, we derive 
the bulk viscosity of the hyperonic matter and discuss 
their role in r-modc damping. In section III, we discuss 
the results of the present calculations. Finally, in section 
IV, we summarize the results and give possible outlook 
of the present investigation. 
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The first two lines of the above Lagrangian density rep- 
resent the interaction of baryons \& b with the aforesaid 
mesons. In the next two lines, we have the kinetic and the 
non-linear terms in the pseudoscalar-isovector pion field 
V, the scalar field 'er', and with x 2 = tt 2 + a 2 . Finally 
in the last two lines, we have the field strength and the 
mass term for the vector field i uj' and the iso-vector field 
'/? meson. The terms in equation (1) with the subscript 
' B' should be interpreted as sum over the states of lowest 
baryonic octet. In this paper, we shall be concerned only 
with the normal non-pion condensed state of matter and, 
so we take < n >= 0. 

The interaction of the scalar and the pseudoscalar 
mesons with the vector boson generates a dynamical mass 
for the vector bosons through spontaneous breaking of 
the chiral symmetry with scalar field getting the vacuum 
expectation value Xo- Then the masses of the baryons, 
the scalar and the vector mesons, are respectively given 
by 



m B = 9aBXQ, 



2Ax , 



lu> — 9uBXQ 



(2) 



II. FORMALISM 



A. Effective chiral model and equation of State 



In this subsection, we briefly summarise the salient fea- 
tures of the chiral model including hyperons, which has 
been considered earlier in the context of neutron star 
matter. It will be used later to calculate the bulk viscos- 
ity. 

The effective chiral model that we shall con- 
sider is a generalisation of the model considered in 
Ref.[5S] to include the lowest lying octet of baryons 
(n,p, A , £~ ,0 ' + , H~'°). They interact via the exchange 
of the pseudo-scalar mesons it, the scalar meson tr, the 
vector meson uj and the iso-vector p— meson. Lagrangian 
density under consideration is given by 



In the above, xq is the vacuum expectation value of 
the a field. We could have taken an interaction of the 
p— meson with the scalar and the pseudoscalar mesons 
similar to the omega meson. However, a dynamical 
mass generation mechanism of the p— meson in a simi- 
lar manner will not generate the correct symmetry en- 
ergy. Therefore, we have taken an explicit mass term for 
the isovector p— meson similar to what was considered in 
earlier works (32] [35] [36] . 

We use mean- field approximation ([23]) to evaluate the 
meson fields in our present calculations. In the mean-field 
treatment, one assumes the mesonic fields to be uniform, 
i.e., without any quantum fluctuations. The details of 
the model that we use in our present investigation and 
its attributes such as the derivation of the equation of mo- 
tion of the meson fields and its equation of state (e & P) 
of the many baryonic system, can be found in our pre- 
ceding work [3B]. The vector meson and the iso-vector 
meson field equations in the mean field framework are 
then given by 
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In the above equations the quantity p B is the baryon 
density and I3B is the third component of the isospin of 
each baryon species. 

The scalar field equation can be written in terms of the 
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where the effective mass of the baryonic species is m* B = 
Ym B and c aB = g 2 aB /m 2 a are the c wB = gl B /m 2 are the 
usual scalar and vector coupling constants respectively. 
Similarly, in the present model describing dense mat- 
ter, the uj— meson mass is generated dynamically. This 
vector meson mass enters in Eq. |5]) through the ra- 
tio Cu; — (g^/m^) 2 = 1/xq 2 . Various parameters of the 
model for the nuclear matter case (i-th meson nucleon 
couplings, C CT , C^, C p and the nonlinear couplings B and 
C) are fitted from nuclear matter saturation properties 
[3"2l 131k 38 . It was shown recently that these parameters 
are rather constrained by the nuclear matter saturation 
properties like the binding energy per nucleon, satura- 
tion density, the nuclear incompressibility, as well as the 
asymmetry energy, the effective mass of the nucleon and 
the pion decay constant f n . Further in Eq.^, the quan- 
tities p B and psb are the baryon density and the scalar 
density for a given baryon species given respectively as, 
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where k B is the fermi momentum of the baryon and 7 = 2 
is the spin degeneracy factor. 

For a given baryon density, the total energy density 
'a' and the pressure 'P' can be written in terms of the 
dimensionless variable Y — x/xq as 
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The terms in eqns. (3) and (4) with the subscript l B' 
should be interpreted as sum over all the states of the 
baryonic octets. The meson field equations for the <r, to 
and p— mesons are then solved self-consistently at a fixed 
baryon density to obtain the respective field strengths. 
The EoS for the j3— equilibrated for the hyperon rich 
matter is obtained with the requirements of conservation 
of total baryon number and charge neutrality condition 
given by [35] 
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where p B and pi are the baryon and the lepton (e,/i) num- 
ber densities with Q B and Qi as their respective electric 
charges. 



B. Hyperon Bulk-viscosity 

Bulk viscosity characterizes the response of the sys- 
tem to an externally oscillating change in the volume. 
The volume expansion or contraction leads to a change 
in density or the chemical potential of the system. This 
drives the system out of chemical equilibrium. The equi- 
librium is restored by the microscopic processes. If the 
equilibrium time scales are comparable to the oscillat- 
ing time scales, there will be energy dissipation. In the 
context of neutron star, the typical oscillation frequen- 
cies are less than a kilo hertz. Therefore, the microscopic 
processes that will be relevant are the weak processes. 

It is already known that the non-leptonic processes 
containing hyperons lead to high values of bulk viscosity 
for rotating neutron stars with temperature 10 9 — 10 10 K 
[TBI fT8l [39] . The leptonic processes are suppressed by 
smaller phase space factors. Thus the relevant reactions 
which are going to give a lower limit on the rates (or 
upper limit on bulk viscosity) are 
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The coefficient of bulk viscosity relates difference be- 
tween the perturbed pressure p and the thermodynamic 
pressure p to the macroscopic expansion of the fluid as 
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where v is the velocity of the fluid element and £ is the 
coefficient of bulk viscosity, which, in general, is complex 
in nature [42] . 

The relativistic expression for the real part of C, which 
corresponds to the damping, can be calculated in terms of 
microscopic equilibrium restoring reaction rates |18l I42j . 
Within a relaxation time approximation, the real part of 
C is given as, 
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where 7oo and 70 are the "infinite" and "zero" frequency 
adiabatic index respectively, uj is the angular frequency 
of the perturbation in co-rotating frame and r is the net 
equilibrium restoring microscopic relaxation time. The 
expression for — 70 is 
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Here ng corresponds to the total baryon density and x = 
— is the neutron fraction. Thus the difference 700 — 7o can 
be calculated from a given equation of state. The angular 
frequency to of the r - mode (1=2, m=2) in a co-rotating 
frame is given in terms of the Keplcrian frequency fix of 



the rotating star as lu = 
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The prominant reactions involving the lightest hyper- 
ons, E - and A , which have higher population in a given 
star are as given by equations (111 and (12). The rates 
of these reactions can be calculated from the tree-level 
Feynman diagrams involving the exchange of a W boson. 

We are not be considering the other reaction n + n < — > 
n + A° since it has no simple W-boson exchange picture. 
We shall discuss more regarding this in Section III. The 
relaxation time r (at a temperature T), when both S _ 
and A are present, is given by [TBI US] 
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where ks is the Boltzmann's constant and k\ and fcs are 
the Fermi momenta of these hyperons. Sfi = 5(i n — <5/iA 
is the chemical potential imbalance. 5x n = x n — x n is 
the small difference between the perturbed and equilib- 
rium values of the neutron fraction. (|.M 2 |) are the angle 
averaged, squared, summed over initial spinors matrix el- 
ements of the reactions calculated from the Feynman di- 
agrams. We refer [T51 [T5] for the expressions for (|.Mf;|) 
and (l-M/J). We note that the contribution from A hy- 
perons will not be present in Equation (17) while con- 



sidering a neutron star medium before the appearance of 
A. 

The factor Sfj,/riB5x n is determined from the con- 
straints imposed by the electric charge neutrality and the 



baryon number conservation given respectively as 

5xp-5x-E = (18) 
Sx n + Sxa + Sx p + 5x-£ — (19) 

together with the condition that the non-leptonic strong 
interaction reaction 



A u 



P 



(20) 



which has a higher rate, is in equilibrium compared to 
weak interaction processes giving rise to the bulk vis- 
cosity. Equilibrium of this reaction implies that both the 
reactions ( 11 ) and ( 12 ) have equal chemical potential im- 
balance, 

Sfi = Sfi n - 5/j,a = 2<5/i„ - 5 n P - 5 fiy,. (21) 
Using these constraints, we can write, 
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These expressions are for the case when both the £~ and 
A hyperons are present. One can not use the reaction 
(20) while considering the region where we have only S~ 



hyperons. In that case, instead of Eq.(22), we have 
2S(i 
nsSxr, 
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Now the ay's can be found out using the expression for 
baryon chemical potential and the equations of motion of 
the mesonic fields given by equations |3p|. In general, 
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As before, Y is the scalar field expectation value in the 
medium in units of its vacuum expectation value. Fur- 
ther, are calculated from the scalar field equation 
Eq.|5]) and are given by 
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with 
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Using equations ( 23p7 ), one can compute the relaxation 
time from Eq.(17l and hence the bulk viscosity given in 



Eq.(15l, for a given equation of state. 



C. R-mode damping 

As mentioned in the introduction section, the r-modes 
correspond to the axial modes where the restoring force is 
the Coriolis force. The r-mode frequency lo c in the coro- 
tating frame, to first order in J7, the rotation frequency 
of the star is given by [30] 
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The r-modes correspond to I = m[4T]. The r-mode fre- 
quency observed by an inertial observer is given by 
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Thus, I — m = 2 mode becomes most unstable to emis- 
sion of gravitational waves due to Chandrasekhar Fried- 
man Shultz (CFS) mechanism. In a rotating star emis- 
sion of these waves causes the modes to grow. This in- 
stability can get damped due to the various viscosities of 
the stellar matter. This happens when the damping time 
scales associated with these viscous processes are compa- 
rable to the gravitational radiation (GR) time scale. 

We need the expressions for the time scales associated 
with the dissipative processes and GR in order to un- 
derstand the nature of the damping of the r-modes. The 
imaginary part of the dissipative time scale (which causes 
the damping) is given by [HI [43] 
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where i labels the various dissipative phenomena like hy- 
peronic bulk viscosity (B), bulk viscosity due to Urea 
processes (U), shear viscosity (rj) GR etc. In the above, 
E is energy of the r-mode in the co-rotating frame. This 
can arise both from velocity perturbation as well as the 
perturbation of the gravitational potential. For a slowly 
rotating star, the dominant contribution is from the ve- 
locity perturbation and is given as 
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with, p(r) being the mass density of the star. Assuming 
the spherical symmetry, mode energy can be reduced into 
an one dimensional integral [5D] as 
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with I — m — 2 for the r-modes and R denotes the radius 
of the star, a is the dimensionless amplitude coefficient of 
the mode, which gets cancelled out in the r calculation. 
This energy is dissipated both by gravitational radiation 
as well as thermodynamic transport of the fluid [T2J H3] • 
The dissipation rate due to the bulk viscosity effects is 
given by 



dE E 
dt 
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Here, in general, |V • Sv\ 2 depends upon the radial and 
the angular co-ordinates. In slowly rotating stars, to the 
lowest order, £ depends only on the radius. Therefore, 
to the lowest order in f2, it is possible to write the bulk 



viscosity dissipation rate in Eq. ( 33 1 as an one dimen 



sional integral by defining a quantity which is the angle 
averaged expansion squared ( |V • Sv\ 2 ). In terms of this 



quantity, Eq. ( 33 1 can be written as 
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where ^|V-<5v| 2 y can be determined numerically 
However Ref.s[lHl US] give an analytic expression 
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which fits to the numerical data. Here G is the gravi- 
tational constant and p is the mean density of the non- 
rotating star. 

Now with the knowledge of density profile p(r) of the 
star, it is straightforward to find out the bulk viscosity 
damping time scales from Equations Eq.(35), Eq.(34|, 



Eq. ( 32 ) and Eq. ( 30 ) , once we know the bulk viscosity co- 



efficient C( r )- ln the case of bulk viscosity time scale aris- 
ing due to hyperons (tb), we can get from Eq.(15) 



. Similarly we can find out the time scale (tjj) associ- 
ated with modified Urea processes, with the help of the 
expression for associated bulk viscosity Qu given by [45] 
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The shear viscosity time scale is given by 
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TABLE I: Parameter set for the model. 



CctJV 


CujN CpN 


B 


C K 


m* N /mN 


(fm 2 ) 


(fm 2 ) (fm 2 ) 


ifm 2 ) 


(fm 4 ) (MeV) 




6.79 


1.99 4.66 


-4.32 


0.165 300 


0.85 



where r\ can be calculated from the prominent nn scat- 
tering and is given by [TO] 



V = 2 X 10 18 ^ 4 r 9 - 2 g/(cm s). 



(38) 



Here p 15 = p/(10 15 g/cm 3 ) and T 9 = T/(10 9 K) are 
density and temperature respectively, casted in dimen- 
sionless forms. Finally the gravitational radiation time 
scale (tgr) is given by [20], 
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The evolution of the r-modc due to dissipative viscous 
effects and GR can be studied by defining the overall 
r-mode time scale T r [TBI |20] , 
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(40) 

It appears in the decay of the mode as e t / Tr and when 
r r > 0, the mode is stable. Now from Equation ( [39] ) 
we can see that tgr < 0, which is indicative of the fact 
that GR allows the modes to grow and drives them to 
instability, while tb, tjj and are positive and thus 
they try to dampen the mode. We can define a critical 
angular velocity He as l/r r (ilc,T) = 0; for a star at a 
given temperature T. Now if the angular velocity of the 
star is greater than Clc then the star is unstable and will 
be subjected to GR emission while stars with angular 
velocities smaller than He will be stable. 



III. RESULTS AND DISCUSSIONS 



Now with the EoS discussed in subsection III Al we set 
out to find the coefficient of bulk viscosity as given by 
equation (151. The parameters in the effective chiral 



model that we have used are given in Table I. The pa- 
rameters were so chosen that they satisfy the constraints 
on the equation of state from the flow data in heavy ion 
collisions [37] . The resulting EOS is plotted in figure [l] . 
To calculate the bulk-viscosity, we first need to calculate 
7oo ^ 70) the difference between fast and slow adiabatic 
indices, from Eq.(16|. This expression can be calculated 
with the help of EoS alone. In figure[2j we plot (700 — 70) 
as a function of the normalized baryon density (n^/no), 
where no = 0.153 fm~ 3 is the nuclear matter saturation 
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FIG. 1: Equation of State used to calculate the stellar config- 
urations 
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FIG. 2: Thermodynamic factor 7^ — 70 that appears in the 
expression for hyperon bulk viscosity, is plotted against the 
normalised baryon density. 



density. The sudden rises in the graph can be attributed 
to the appearance of hyperons with increase of baryon 
density at the cost of neutron number density n n . 

Since we have S _ and A hyperons formed in the system 
with lowest threshold densities, we consider the non lep- 
tonic reactions represented by the Eq. (Ill and Eq.(12) 



and calculate the relaxation time as given by the Eq. 
(17). We note that we have not considered the reaction 



Eq.(13) and further there are several reactions which are 
going to contribute to the net reaction rate [EH [21] • The 
reason for this is that the rate for the process given in 



Eq.(13) is estimated to be an order of magnitude higher 
than the process given by Eq.(12). This leads to sub- 



dominant contribution to the relaxation time and hence 
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FIG. 3: Relaxation time r (in seconds) for the non-leptonic 
processes causing hyperon bulk viscosity is plotted for various 
temperatures T 



FIG. 4: Hyperon bulk viscosity £ in units of g/ (cm s) is plot- 
ted as a function of normalized baryon density for various 
temperatures 



to the bulk viscosity [35]. Thus what we are calculating 
is a lower limit of the net rate which will correspond to 
an upper limit on the bulk viscosity. The matrix ele- 
ments are calculated with the values of Fermi momenta 
and effective masses of various baryonic species obtained 
from the EoS. Here we use axial- vector coupling values 
fjnp — —1-27, (fyA = —0.72 and g n s — 0.34 measured in 
/3-decay of baryons at rest and Fermi coupling constant 
Gp = 1.166 x 10~ n MeV" 2 and sin6» c = 0.222 (where 
6c is the Cabibbo weak mixing angle) (46]. Then, we cal- 



culate 



from Eq. ( 22 ) for the densities where both 



the hyperons are present (n^/no > 2.36), and, from Eq. 
( 23 1 for lower densities where there is only S-hyperon 
present (n B /n = 1.86- 2.36). Further, Eq.p4|-Eq.p7| 



are evaluated using the EoS under consideration. We can 
thus calculate the relaxation time for relevant tempera- 
tures. We show the calculated behaviour of relaxation 
time (in seconds) with temperature in Fig. [3] It is clear 
that the relaxation time increases considerably with the 
decrease of temperature. For a given temperature, the re- 
laxation time is seen to decrease when both the hyperons 
are present, as compared to the case of presence of a singe 
species of hyperons (£), since in this case r value will be 
less according to Eq.(17). We then compute the coeffi- 



cient of bulk viscosity responsible for the mode damping 
in neutron stars from the expression (15). The value 



of maximum frequency is the Keplerian angular fre- 
quency of the rotating star and is set by the onset of mass 
shedding from the equator of the star. The bulk viscos- 
ity coefficient is calculated for the relevant temperatures 
and is plotted against the normalized baryon density in 
Fig. |4j The behaviour of the hyperon bulk viscosity is 
similar to that of the corresponding relaxation time as is 
expected from Eq.(15). The high value of the bulk vis- 
cosity coefficient at the temperature 10 9 K is indicative 




FIG. 5: Density of the star, p (in units of g/cm 3 ), as a function 
of distance r (in km) from the centre (r= 0) to the radius of 
the star (r = R). Threshold densities corresponding to the 
formation of E~ and A hyperons are also plotted. 



of the fact that hyperon bulk viscosity plays a major role 
in the suppression of the r-modes. We note that our bulk 
viscosity values are order of magnitude less than the val- 
ues obtained by [18]. It could be due to the fact that 
unlike their work we are not considering the effect of hy- 
peron superfluidity in this calculation. It might also be 
noted that the non-superfhiid hyperonic bulk viscosity 
calculated in [48] uses an EoS based on a model, where 
only A hyperons are present at the relevant density. 

We next study the effect of hyperon bulk viscosity 
on the r-modes. Here we need to calculate the dissipa- 
tion time scale due to hyperon bulk viscosity as well as 



I 1 1 1 1 1 > 1 

10 7 10 8 10 9 10'" 10" 10 12 

T(K) 

FIG. 6: The temperature dependence of damping time scales 
(in seconds) due to hyperonic bulk viscosity tj, modified Urea 
bulk viscosity tb and shear viscosity t v . tgr represents the 
temperature independent gravitational radiation time scale. 
(The star is considered to be rotating with the Kepler fre- 
quency Qk here). 



0.02 1 — — — 

10 8 1 9 10 10 1 0" 

T(K) 

FIG. 7: Critical angular velocities (normalized to the Kepler 
frequency Qk = 3998 Hz) for a neutron star with mass of 
1.66 Mq is shown as a function of hyperon core temperature. 
The shaded region represents the majority of the observed 
LMXBs. 



due to other dissipative phenomena. If this time scale 
is greater than the gravitational radiation time scale, 
then the r mode is not stable and suppressed. In or- 
der to calculate the dissipative timescales from Equations 
( 30 )-( 39 1 we need to know the density profile p(r), of 
the neutron star under consideration. We need to know 
the Kepler frequency of the rotating star also. We use 
Tolman-Oppenheimer-Volkoff equations to construct the 
non-rotating stellar configurations. The maximum mass 
of the neutron star in this case is found to be 1.65 Mq 
with a radius of 16.7 km. We use Hartle's slow rotation 
approximation to calculate the global properties of ro- 
tating neutron star|49j. We get the maximum mass and 
radius (i?) of the rotating star to be 1.66 M and 18.9 
km respectively. The Kepler frequency in this case is 
found to be £Ik — 3998 Hz. A typical density profile of 
the rotating star is shown in the Fig. [5] . This profile cor- 
responds to a central density of 7.48 x 10 14 gm/cm 3 . We 
have also indicated the densities corresponding to the ap- 
pearance of the hyperons, i.e., threshold densities of both 
E and A hyperons in the graph. From centre of the star 
upto a density of p — 6.34 x 10 14 gm/cm 3 , we have the 
presence of both the hyperons in the star (i.e., up to a dis- 
tance 2.5 km from the centre). The presence of £ alone 
is there upto p = 5.1 x 10 14 gm/cm 3 (another 1.7 km) 
making a hyperon core of radius 4.2 km in the neutron 
star. Hyperon bulk viscosity time scale, and, hence its 
effects on r-mode is very sensitive to the hyperonic core's 
constituent structure and its radius. 

For the rotating neutron star with a mass of 1.66 Mq 
and Q.k = 3998 Hz as considered above, we next eval- 
uate the various dissipative time scales associated with 



the r-mode damping. The dissipative time scale of hyper- 
onic bulk viscosity, denoted by tb can be calculated from 
equations ( 30 >-( 35 ) with the help of density profile of the 
star. Here hyperonic bulk viscosity (£) as a function of 
radius is obtained from our previous calculations of £ for 
the EoS together with the knowledge of stellar density 
profile i.e., £(p(r)). The time scale associated with modi- 
fied Urea processes tu, is calculated in the same manner 
as for tb , by using the equation ( [36| instead of hyperonic 
bulk viscosity. Next, we estimate the shear viscosity dis- 
sipative time scale t v using the equations ( 37 1 and ( 38 1 . 
Finally, the gravitational radiation time scale tqr asso- 
ciated with the r-mode can be calculated with the help 
of density profile using equation (39). Figure [6] shows the 
calculated time scales as functions of temperature, for the 
star rotating with Q.k- From figure [6j we observe that 
in the non-supcrfiuid hyperonic matter, r-modes get sub- 
stantially damped due to hyperonic bulk viscosity only 
at low temperatures (T < 10 s K), whereas the modified 
Urea bulk viscosity suppresses the r-modes rapidly only 
at high temperatures T > 5 x 10 10 K. The role of shear 
viscosity in suppressing the modes is not prominent in 
this temperature range. Consequently, the effect of r- 
mode instability will be prominent in the temperature 
window (10 8 — 5 x 10 10 ) K. The hyperon bulk viscosity 
suppresses the instability for temperatures below 10 8 K 
while modified Urea processes suppress the instability be- 
yond 10 10 A. 

Now we are in a position to calculate the critical an- 
gular velocity VLc of the neutron star. f2c is obtained 
by solving equation (40); - ^ c T ^ = 0, for a particular 
value of T. At this frequency, the energy fed into the 
r-mode per unit time by gravitational radiation is equal 
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to the energy dissipated per unit time. A star rotating 
above this critical frequency will be subjected to r-mode 
instability. We have shown the fi^ (scaled to fix) in Fig. 
[7] for the temperature regime of interest. Since Qk de- 
termines the maximum allowed rotation rate for the star, 
stable rotation at any temperature will have Qc I = 1 
as an upper bound. In this figure, the region above flc 
curve is unstable and a star rotating in this region will 
be rapidly spun down to an angular frequency below tie- 
As expected the instability window exists in the temper- 
ature regime (10 8 — 5 x 10 10 ) K, where gravitational radi- 
ation is dominant and not suppressed substantial, which 
shows that the neutron star with hyperonic core is un- 
stable in this region. Low temperature regime hyperonic 
bulk viscosity damps the r-mode effectively whereas nu- 
cleon dominated modified Urea bulk viscosity is the cause 
of mode damping at high temperatures. The minima of 
flc curve occurs at T sa 5 x 10 10 K with Q c ~ -04 Ojf, 
which is indicative of the fact that the r-mode instabil- 
ity is rather strong in this hyperon core scenario. The 
shaded box in the Fig. [7] is where most of the observed 
Low Mass X-ray Binaries (LMXBs) are found. They have 
a core temperature in the range of (2 x 10 7 — 3 x 10 8 ) K 
with rotation rate between 300 to 700 Hz (TTJl [5U]. In 
our case, LMXBs are placed in the stable region, unlike 
conventional neutron stars with npe matter |43) . 

IV. SUMMARY 

Rotating equilibrium configurations of self gravitating 
fluids are subjected to various possibilities of instabili- 
ties at large rotation periods. In the present work, we 
have investigated the r-mode instability which is known 
to limit the angular velocities of rapidly rotating stars. 
The r-mode and the related instabilities are damped by 
various viscosities of the matter in the interior of the 
neutron star. Thus the microscopic models describing 
the matter in the interior of the star get constrained by 
the observations of the rapidly rotating pulsars. 

In the present work, we have confined our attention to 
the case of neutron star with a hyperonic core. For the 
description of the matter in the core of the neutron star, 
we have used an effective chiral hadronic model gener- 
alised to include the lowest lying octet of baryons. The 



parameters of the model are chosen that are consistent 
with the flow data in heavy ion collisions, nuclear matter 
properties as well as observation of high masss neutron 
stars. In the present work, we have computed the coef- 
ficient of bulk viscosity due to the the hyperonic matter 
in the core of a neutron star and the resulting effects 
on the r-mode instability. It turns out that hyperon 
bulk viscosity within the model is effective in damping 
the instability for temperatures below 10 8 K . Beyond a 
temperature of about 10 10 K, the bulk viscosity due to 
modified Urea processes become effective in damping the 
r-mode instability. Shear viscosity of hadronic matter 
becomes effective in damping only at low temperatures. 
Within the model it turns out that the the bulk viscosity 
in normal hyperonic matter does play an important role 
in spinning down fast rotating neutron stars. However, 
superfluid hyperonic matter or quark matter in the core 
can change this conclusion. 

We have not considered in the present work the phase 
transition to quark matter which could be most likely in 
a color superconducting phase. The role of quark matter 
in the context of r-mode characteristics has been dealt 
with in Ref.[in] which shows that the r-mode instabil- 
ity gets suppressed by both normal quark matter as well 
as gapped quark matter in the color flavor locked phase. 
Further, the neutron stars are also endowed with strong 
magnetic fields the effects of which on the r-modes have 
not been included in the present work. Ultra strong mag- 
netic field seem to increase the instability window for 
quark matter It will thus be interesting to examine 
the scenario with a phase transition to quark matter and 
the effect of magnetic fields in the context of r-modes 
for a hybrid star with a crust of hadronic matter. Work 
in this direction is under progress and will be reported 
elsewhere |51j . 
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